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Abstract 

We derive a general formalism for evaluating the high-frequency limit of the thermoelectric 
power of strongly correlated materials, which can be straightforwardly implemented in available 
first principles LDA+DMFT programs. We explore this formalism using model Hamiltonians 
and we investigate the validity of approximating the static thermoelectric power So, by its high- 
temperature limit, S*. We point out that the behaviors of S* and So are qualitatively different for 
a correlated Fermi liquid near the Mott transition, when the temperature is in the coherent regime. 

O 

When the temperature is well above the coherent regime, e.g., when the transport is dominated 
by incoherent excitations, S* provides a good estimation of So- 
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I. INTRODUCTION 



Thermoelectric energy harvesting, i.e. the transformation of waste heat into usable elec- 
tricity, is of great current interest. The main obstacle is the low efficiency of materials for 
converting heat to electricity [HE]- Over the past decade, there has been a renewed interest 
on thermoelectric materials, mainly driven by experimental results [3]. 

Computing the thermoelectric power (TEP) in correlated systems is a highly non-trivial 
task and several approximation schemes have been used to this intent. The well-known 
Mott-Heikes formulajl] gives an estimate of the high temperature limit of TEP [3] in the 
strongly correlated regime. A generalized Boltzmann approach including vertex corrections 
has been developed in Ref. [6] and applied to several materials. Thermoelectric transport 
at intermediate temperature was carefully investigated in the context of single-band and 
degenerate Hubbard Hamiltonians, by dynamical mean field theory (DMFT) [7J [8] . Kelvin 
formula was also revisited for various correlated models in Ref. [9] very recently. 

The high frequency (AC) limit provides another interesting insights to gain further un- 
derstanding of the thermoelectric transport in correlated materials, and is the main interest 
of this work. The thermopower in the high frequency limit of a degenerate Hubbard model 
near half-filling was considered in Ref. [8j, where the authors generalize the thermoelectric 
response to finite frequencies in the high temperature limit. The same limit was studied 
recently by Shastry and collaborators, who have developed a formalism for evaluating the 
AC limit of thermoelectric response using high temperature series expansion and exact di- 
agonalization. The methodology was applied to a single band t-J model on a triangular 
lattice [TU1 ITT] . The authors pointed out that the AC limit of TEP (S*) is simple enough 
that it can be obtained by theoretical calculations with significantly less effort, while still 
provides nontrivial informations of the thermoelectric properties, and give an estimation of 
the trend of So. 

In this work, we investigate the high frequency limit of TEP, S*, by deriving an exact 
formalism in the context of a general multi-band model with local interactions. We show 
that S* is determined by the bare band structure and the single-particle spectral functions. 
The relation between the conventional TEP, i.e., obtained at zero frequency (So) and the AC 
limit S* is discussed from general arguments on the single particle properties of correlated 
systems at low and high temperatures. The analytical derivation of S* is compared with the 
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frequency dependent thermopower of the one band Hubbard model, solved by dynamical 
mean field theory (DMFT) on the square and triangular lattices. The formalism derived 
in this work can be conveniently implemented into first-principles calculations of realistic 
materials, such as in the LDA+DMFT framework [T2l [T3] . 

This paper is organized as follows. In Sec. II A general formalism of dynamical ther- 
moelectric transport coefficients is summarized to define the notation. In Sec. II B[ exact 
formulae to evaluate S* are derived for a general tight-binding model with local interactions. 



In Sec. Ill, we apply the formalism to one-band Hubbard model on square and triangular 
lattice. The low and high temperature limit behaviors of S* are discussed and compared to 



those of Sq. Numerical results are presented in Sec. IV Sec. M summarizes the paper. 



II. DYNAMICAL THERMOELECTRIC TRANSPORT FUNCTIONS AND HIGH- 
FREQUENCY LIMIT OF THERMOPOWER 

A. General formalism 

Electrical current can be induced by gradient of electrical potential and temperature. 
The phenomenological equations for static(DC limit) external fields are [14] 

Jf = (-^) + L n (v4) > C 1 ) 

J2 = LTi (-^V„/i) + L'22 (V^) • (2) 

We only consider the longitudinal case. and are x ~ component of particle and 
heat current, respectively. V x jl and V ' x \ are generalized forces driving Jf and Jf. jl = 
— eV, in which fi is chemical potential and V is the electric potential. Lf? are transport 
coefficients. We follow the definition in Ref. [H], which explicitly respects the Onsager 
relation, Lf? = L^f. Transport properties can be defined in terms of L??. For example, the 
electric conductivity a, thermoelectric power S, and the thermal conductivity k are 

° = fUn, (3) 

1 T XX 

S = -4t^, (4) 

, r, (-^12) 2 A /c\ 

^22 - —F^x- ■ ( 5 ) 



y AA J^xx 
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In following context, we use ks — e — h — 1. The practical value of S is recovered by 
multiplying the factor fce/e = 86.3fiV/K, which we use as unit for thermopower. 

In conventional thermoelectric problems, L l J x is theoretically defined and experimentally 
measured at the DC limit. The extension to dynamical (frequency) case is absent in standard 
textbooks but has been studied in detail in Ref. [10J. Here we give the outlines of the for- 
malism. Borrowed from Luttinger's derivation [15], an auxiliary "gravitational" field coupled 
to energy density is defined. An "equivalence" between the fictitious gravitational field and 
the temperature gradient is proved. Then the transport coefficients L?? can be written in 
terms of correlation functions between particle current and(or) energy current. In Ref. [10], 
this formalism is generalized to temporally and spatially periodic external fields, thus the 
transport coefficients become momentum- and frequency-dependent functions, L? x (q, u). 

Some interesting remarks can be made on Lfj E (q, ui). The thermodynamic limit corre- 
sponds to q — > 0, and the static fields correspond to the tu — > limit. The L xx in Eq. (|T| and 
Eq. ([2]) can be approached by the fast limit, i.e., taking q — > first and then u — > 0. If we 
define the "phase velocity" of the external field, v = ^ , the fast limit means v — » 00, which 
gives the name "fast". The slow limit means reserving this order, u> — > first and then q — >, 
thus v — > 0. The slow limit gives the Kelvin formula of thermopower discussed in Ref. [9]. 
The high-frequency(AC) limit means u> — > 00. In this case, we take the thermodynamic 
limit, q — > first, and then 00 — > 00. But from the general formalism in Ref. [10], it can be 
shown that the order of taking limits does not matter. 

The dynamical transport coefficients with q — > are given by, 

poo p/3 

Ll x (u)=T dte^ M0+)t dr{J*{-t-iT)J?). (6) 
Jo Jo 

For a given Hamiltonian H, the current operators are defined by following the conserva- 
tion laws [14]. 

dO x 

Jr = -Q±- = i[H,Of\. (7) 
Of is the x-component of particle and heat polarization operator. Specifically, 

Of = J>?n,, (8) 

i 

O x 2 = J2 R i - A*»i) , (9) 

i 

where rij and hi are local particle and energy density operators. The explicit forms of rij and 



hi are determined by the Hamiltonian of specific models. In next subsection, we will write 
Oi and give Jj for a general multiband model. 

At DC limit, the imaginary part of L x ^{uj = 0) is zero, thus So is determined by the real 
parts. For convenience, define 

L% = ReL^(O), (10) 



then we have 



1 L 



o 



So = ReS(u = 0) = --^£. (11) 
At AC limit, L^(u) is dominated by the imaginary part, with a 0(1/ uS) leading order, 



\™L*J{u) = -L* + O{- 2 ). (12) 

Using Lehnman's representation, it has been shown that L|- defined above is, up to a factor 
of i, the expectation values of commutators between current and polarization operators jH 
EUim, i-e., 

l%=i([j;,Of\). (13) 

Consequently, TEP at AC limit is 

S * = ReS (u oo) = (14) 

L*j can be related to ReLjj(cj). Applying Kramers-Kronig relation and keeping the 
leading order in 1/ou, we have 

1 f°° 

L*j = / cL,ReL™(uj). (15) 

Thus L*j is also connected to the sum rules of dynamical quantities. For example, L* x is 
proportional to the sum rule of conductivity [161 E]- 

L* = - / dwRea(w). (16) 
^ Jo 

Other sum rules are also derived in Ref . [10] and 1 1 1 1 . 



B. General formula of L*, 



Now we explicitly evaluate the commutator in Eq. (13) for a general tight-binding Hamil 



tonian with local interaction, which will determine the AC limit of TEP in this system. We 
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start with the following Hamiltonian 

# = - E C c L c ^ + E e ^ c l c 



t/i 



(17) 



i, j are site indices, a, /3, /x and z/ denote local orbitals. is the hopping integral, and 
Uaf3fiu is the matrix element for Coulomb interaction between local orbitals. is energy 
level of local orbitals. The particle polarization operator is 



^ = E^E c L c ^ 

i fi 

and the heat polarization operator is 



0\ = 



The current operators turn out to be 

Jt = i[H,Of\ 



and 



J 2 * = i[if,0*] 



-\ E - ^) + e ^ - 2 ^ 



E - 



X 



a'fi'u' 
a'P'v 1 



(18) 



(19) 



(20) 



(21) 
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In the literature [18J, Jf is also written in a more compact form using the equation of 
motion in Heisenberg picture, 



^2 — 2 fai Ri ) tfj \fin c i v °l^jv j ■ 



in which the dot means the time derivative, 

To compute L\ x and L* 2 , we need to further evaluate the commutators between current 
operators and polarization operators. For L* x , this is simple and straightforward, 

However, L\ 2 leads to a complicate formula, 



+^ E C - R i) («* + e »- ^ <i c ;*> 

ij,fiv 

+\ E - R i) 

ij,fj,v 

X I {Uva'n'v' — U a i u ^ u i) (cJ M ct ,CjvCjy) 

\a'fj,'u' 



(23) 



But this formula can be significantly simplified if we look at the equation of motion for the 
following Greens's function, 

GJ{t) = -(T T c jv {r)cl). (24) 

T T is the time-ordering operator in imaginary time. Its equation of motion reads, 

dGT(r) 



— {U a i vl ji v i — U ua 'fj,i u i) 

X ( T r C Ja' ( T ) C JV' ( T ) C jV ( T ) C l) • 
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Taking the r — > limit leads to 



y ^ji \ T ) ST^ ,VV> I t \ 



/^)(4^). (25) 



Substituting the last term in Eq. (23) by the right hand side of Eq. (25), we get 
1 



'\2 



£ W ( R * ~ R i) 



(R? - R*f - (Rj - Rf) 2 ] (clc jv ) 



EC(^-^) 2 Hm |W(t). (26) 



o- dr 



Using the fact that 



and performing Fourier transformation in both real space and imaginary time, we get 

L *n = -o £^°~ E V-W-) G ^n\ (27) 

u„ k,/iu ^ x ' 

and, 



'VI 



P 2 ^ u ^ V 9k J \ dk x 



'9e^\ . /<9 2 ef 



Gr(^n)- (28) 



is Fourier transformation of hopping amplitudes, 



fll/ 



J2e ikR t» u (R), 



(29) 



where we have utilized the translational invariance, 



(30) 



It is straightforward to convert the Matsubara summation to the integration in real 
frequencies. 



l* u = r dwj^i 

J -oo ,,,, V 



dk 2 

k.fiu v x 



(31) 
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and 



L 



12 



E 



dk x 



dkl 



(32) 



/(w) = 1/(1 + exp(/3a;)) is the Fermi function. A^iu)) = — -G v ^{u}) is the spectral function 



Eq. (27), Eq. (28), Eq. (31) and Eq. (32) are main results in this work. They are derived 



from a general formalism of dynamical thermoelectric transport outline in Sec. |II A| and a 



multiband Hamiltnian, Eq. (17). The equation of motion is exact and no approximation 



is assumed in the derivation. These equations indicate that L\ x and L* 2 , and thus S* are 
determined by the non-interacting band structure and the single-particle spectral fundtion. 



III. S AND S* IN A ONE-BAND HUBBARD MODEL 

In this section, we discuss Sn and S* of one-band Hubbard model in the scenario of 



dynamical mean field theory (DMFT), using the formalism we presented in previous sections. 
The Hamiltonian of one-band Hubbard model is 

ij,a i 

In DMFT, it is mapped to a single-impurity Anderson model [TS] supplemented by the self- 
consistent condition, which reads, 

1 



52G k (zu> n ). (34) 



On the left hand side is the local Green's function on the impurity. A(iu n ) is the hybridiza- 
tion function of the impurity model. On the right hand side, Gk{ioj n ) is the Green's function 
of lattice electrons, 

G k (iuj n ) 



iu n + /i - e k - £(zu; n ) 

with tfc the non- interacting dispersion relation of the lattice model, and Ti(ioj n ) the self energy 
for both local and lattice Green's function in the self-consistent condition. In DMFT, both 
coherent and incoherent excitations in a correlated metal are treated on the same footing[20j. 

In DMFT, the evaluation of transport coefficients,e.g., Eq. (|6]), can be significantly sim- 
plified. Because the fc-dependence falls solely on the non-interacting dispersion e k , the vertex 



9 



corrections vanishes|21j. Consequently, ReLy(o;) can be written in terms of single-particle 
spectral function in real frequency. 

2 



H,eLij(uj) 



dk x 



ff(u')-f(cu' + cu)\ 
x f — ^ '- \ A k (w')A k (u/ + u). 

(35) 

Notice that here the dependence of HeLij(u) on the single-particle spectral function is 
generally approximate for a finite-dimensional system, which is achieved due to the vanishing 
of vertex corrections exact only in infinite dimensions. But the dependence of L*- on single- 



particle spectral function is exact, as pointed out at the end of Sec. II A 



Another question is on the sum rule of the approximate HeLij(u), i.e., if we substitute 



Eq. (35) into the definition of L*-, Eq. (15), wether or not it will give the same form of L*- as 



we have derived in last section. The answer to this question is yes and we a brief proof for this 
one-band case in the Appendix but the extension to multiband case is straightforward. This 
means that ignoring vertex correction will modify the distribution of weight in ReLjj(w), 
but will not change the integrated weight. 

The DC limit of HeLij(u), L°- can be obtained by takeing the limit u — > 0, which gives, 



de k 

dkrr 



dojoj i+j - 2 



dui 



AM 



(36) 



Therefore in the framework of DMFT, So is computed from Eq. (36). The AC limit, 



S* can be computed from Eq. (27) and Eq. (28), or Eq. (31) and Eq. (31). In principle, 



Matsubara frequency and integration over real frequency give identical results. But in 
practice, especially in numerical computations on correlated systems, correlation functions 
in Matsubara frequencies are more easily accessible. For example, among various impurity 
solvers in DMFT, quantum Monte Carlo method(QMC), i.e., Hirsch-Fye method|22j and 
recently developed continuous time QMC[23, 24] are implemented in imaginary time. To 
get correlation functions in real frequencies, numerical realization of analytical continuation 
has to be employed, such as maximum entropy method, which is a involved procedure and 
usually special care has to be taken of. In this case, a formulae in Matsubara frequencies 
will significantly simplify the calculation. 
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Due to the bad convergence of the series, Eq. (27) and Eq. (28) are not appropriate for 
direct implementation into numerical computations. Following standard recipe(separating 
and analytically evaluating the badly convergent part), we transform them in a form more 
friendly to numerics. For the one-band Hubbard model, 

*;.=e(^) ft E -i). < 37 > 



and 



L\ 2 = S(fjr) lj2 ReGk ^ x [l + 2w„ImG fc (ia; n )]. (38) 

A. Low temperature limit. 

At low temperatures(low-T), the derivative of Fermi function, (—df(uS)/duS) in the inte- 



grand of Eq. (36) becomes Dirac-5 function- like, thus only the low energy part of the spectral 
weight near Fermi surface contributes to the integral. The low energy part of the self energy 
of a Fermi liquid £(a>) can be approximated by a Taylor expansion in terms of u and T. 



Re^U-} ~ ( 1- — 1 u.' 



Z z Z A 

Previous studies[3 showed that at low-T limit, L U oc Z 2 /T and L\ 2 oc ZT, thus S = 
-L° 12 /(TL° n ) oc T/Z. 

Since we are interested in the relation between So ad S*, it would be convenient to 
write L* l2 an d i n terms of the conventional transport function, (dek/dk x ) 2 . This can be 



achieved by performing integration by part on the summation over k in in Eq. (31) and Eq 



32 , then we have 

L ij ~ L ij,i + L ij,in 



11 



with 



df(u) 
dcu 



xu^- 2 [--)l m [G k (uj)Z(tu)}, (39) 



7T 

2 



xlm 



G( £)U )|-(^- 2 (l-Z( w ))) 



(40) 



where we have defined 



Z(u) 1 



l-dE(u)/du 

We introduced the function Z(oj), which is dependent on the derivative of self energy 



with respect to energy oj. The integrand in L* j 7 (Eq.(39)) also has the derivative of Fermi 
function. Also notice that at low-T, Z(u = 0) = Z, which is the renormalization factor of 
correlated Fermi liquid. Then L*-j resembles L\ 2 except for the power of IraG^w). Low 
temperature expansion show that L\ XI oc Z, and L* 2I oc T 2 . Therefore, if L* llH and L* 12II 
were absent, S* = — {TL\ 21 )jL\ ll oc T/Z, which is similar to the low-T behavior of Sq. 
However, L* 1H and L* l2II do not vanish in general at low-T limit. First, at low-T limit, 



the integral over u> in Eq. (40) 



/oo rO 
dujf{uj) is replaced by / du. 
■oo J — oo 



Then both the real and imaginary part of Gk{uj) and Z(oS) below Fermi surface have to 
contribute to the leading order of L*j n , unless is independent, or at least weakly 

dependent on u, leading Z(u) ~ 1, and then the integrand in L*- n would vanish. But this 
in general can not be true. For example, in a correlated Fermi liquid phase near the Mott 
transition of Hubbard model, E(w) contains the information of coherent quasiparticles at 
Fermi surface as well as that of incoherent excitations in high-energy Hubbard bands, thus 
S(o;) will depend on u in very different ways at these separated energy scales. At low energy 
scale, Z(oj ~ 0) ~ Z, and Z is significantly less than 1 near Mott transition. Therefore, at 
low-T limit, L*- n will exhibit a finite value at low-T limit. So the total value of L\ 2 will be 
dominated by L* l2II instead of the ~ T 2 contribution from L\ 2I . The finiteness of L\ x can 
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be also justified by the general sum rule Eq. (31), which indicates that L\\ is proportional 
to the kinetic energy. Consequently, S* will diverge 1/T-like at low-T limit for a correlated 
Fermi liquid. 

There are some circumstances in which Z(u) = 1 and L\^ u vanishes. One example is 
that in a static mean field theory, such as Hartree-Fock approximation, £(o>) is independent 
on u), thus in static mean field theory, it is possible that S* can show a similar behavior to 
that of Sq at low temperature. 



B. High temperature limit. 

In the literature, the high temperature limit of thermopowerjl], or known as Mott-Heikes 
formulor, has been widely used as a benchmark for thermoelectric capability [5 J for correlated 
materials. Here we discuss the high temperature limit of S* implied from the formulae we 
have derived. 

The high temperature limit relevant for correlated systems was approached by first taking 
the limit U — > oo, which excludes the double- occupancy in hole-doped systems or the 
vacancy in electron-doped systems, then taking the high temperature limit T — > 0. This 
leads to two major simplification. First, by definition in thermodynamics, 

a = -( ds \ 

T \dN) EV 

Here s is the entropy and N is number of electrons, s can be calculated by counting all 
possible occupation states satisfying the U — Y oo limit. It turns out that ^ is a constant 
determined by the electron density. Thus /i is proportional to T at high temperature. The 
second simplification is that at high temperature, we can approximate the single particle 
spectral function by a rigid band picture, namely, 

A k (u) = A k (u-n). (41) 

Ak(u)) is a function of uj but independent of temperature and chemical potential. Applying 
these simplification to Eq. (31) and Eq. (32), and keeping the leading order in T, we have 
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Therefore, at high temperature limit, 

S* = -3^ = ^, (42) 

This is the same result to the high temperature limit of So in Ref. [3]. Thus the leading 
order of S* is identical to the leading order of S at high temperature. 



IV. NUMERICAL RESULTS 



In this section, we compute the dynamical thermoelectric power S(u) by dynamical 
mean field theory (DMFT). We use exact diagonalization(ED) as the impurity solver. The 
advantage of the ED solver is the Green's functions can be computed simultaneously in real 
and Matsubara frequencies. Thus we have two approaches to compute the AC limit S*. 



The first one is to substitute the Green's function in Matsubara frequencies into Eq. (27) 



and Eq. (28). The second method starts from computing ReLn(w) and ReLi2(w) from 
spectral functions A/,(u) using Eq. (35) for a wide range of u, Kramers-Kronig relation is 
implemented to compute Imlu(w) and Imi 12 (w), and finally with the value of Z/J4 and L\ 2 



obtained by fitting Eq. (12) at the to — > oo limit. The second method is more laborious but 
here we use it as a check for our formulae in Matsubara frequencies. 

We study one-band Hubbard model on square and triangular lattices and consider only 
the hopping between nearest neighboring sites. 



A. Square lattice 

In this section, we compute the thermoelectric transport coefficients and thermoelectric 
power for a hole-doped Hubbard model on square lattice. We use the bandwidth D as the 
unit for frequency u, temperature T and interaction strength U. For square lattice, D = 8\t\, 
t is the hopping constant. 

In Fig. [T] we show the frequency-dependent quantities for U = 1.75D and n = 0.85. 
Fig. [l]-(a) and -(b) show the thermoelectric transport coefficients Ln(w) and Lvziw) by their 



real(red line) and imaginary part (black line). The real parts are computed from Eq. (35). 
The imaginary parts are computed from Kramers-Kronig relation. Three contributions are 
recognizable in ReL n : i) The low frequency peak due to transition within the resonance 
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peak of quasiparticles. ii) The transition between quasiparticles and the lower Hubbard 
band, which accounts for the hump at w ~ 0.5D. iii) The weight around u ~ U, which 
is due to the incoherent excitations between Hubbard bands. Same features also exist 
in ReLi2(w), but the feature near u ~ i.e., transition between quasiparticles, and the 
transition between quasiparticles and lower Hubbard band, are much less obvious. This is 
because the DC limit L° 2 is dominated by the particle-hole asymmetry of the band velocity 
dek/dk x and the spectral function A^lo), due to the uj l+ i~ 2 = u term in the integrand of 



Eq. (36) for . Thus at small u , ReL 12 (o;) is significantly impaired, compared to ReL n (w). 
Therefore the transition by incoherent excitations around u ~ U takes a major part in the 
total weight in ReLi2(w), and the sum rule of HeL^co) , i.e., L\ 2 , is also dominated by the 
incoherent excitations. Imln(u) and ImLi 2 (u;) are odd functions of u and vanish at u = 0. 
It is evident that the real parts approach to zero much faster than the imaginary parts at 
AC limit(w — > oo),. Fig. [l]-(c) shows the evolution of ReL 12 (w) as temperatures. The 
dominance of the incoherent excitations is robust as the variation of temperature. Fig. [l]-(d) 
shows the real part of thermoelectric power, ReS'(a;) for T = 0.0625D and T = 0.0875-D. 
The inset blows up the region near oj = 0, indicating that So displays + or — signs at 
different temperatures. 

In Fig. [2]- (a) and (b) we show So and S* at various temperatures. On the one side, in 
Fig. [2]- (a), 5*o presents multiple changes of sign with temperature increased. The sign change 
at lower temperature(T ~ 0.1D) demonstrates the crossover from the low-temperature hole- 
like coherent quasiparticles to incoherent excitations at intermediate temperature. Around 
T = 0.2D, where So reaches its maximum positive value, where the coherent quasiparticles 
have almost diminished. The second sign change around T = 0.6D indicates a subtle 
competition between the spectral weight of lower and higher Hubbard band. As temperature 
increases, the asymmetry between the two Hubbard bands near Fermi surface becomes less 
significant because more spectral weight from the higher Hubbard band takes part into the 
transport and the sign of 5o is determined by the difference between the weight of lower 
and higher Hubbard. This crossover is thus considered to be responsible for the second sign 
change[8] and also has been observed experimentally [26J. Therefore, above T = 0.6D, the 
transport is completely dominated by incoherent excitations from both Hubbard bands. On 
the other side, in Fig. [2j-(b) , the situation for S* is quite different. S* does not change sign 
and keeps negative in the shown temperature range. Towards low temperature, S* blows up, 
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FIG. 1: Frequency-dependent transport coefficients and thermoelectric power of a hole-doped one- 
band Hubbard model on square lattice. U = 1.75D and n = 0.85. (a) ReLn(o;) and ImLn(w) at 
T = 0.1251?. (b) ReLi 2 (w) and ImL 12 (w) at T = 0.125D. (c) The evolution of KeL 12 (uj) with 
temperature, (d) Re<5(£j) at T = 0.0625Z? and T = 0.0875D. The inset blows up the region near 
oj = 0. 



consistent with our argument based on a Fermi liquid self energy in Sec. HI A Towards high 
temperature, i.e., when the temperature is well above the coherence regime, So and S* have 
the same sign and similar magnitude. We notice that Sq in Fig. |2]-(a) does not converge to 
the value predicted by the Mott-Heikes formula in the correlated regime(SMH — l-04k B /e, 
from Eq. (11) in Ref. [1]). This is because in our case, with U = 1.5D, the requirement for 
\t\ T U can not be satisfied for a wide range of temperature. Thus at high temperature, 
e.g., when T > 0.6D, the states with double occupancy can not be excluded and they are 
responsible for the second sign change in So as discussed above. 

In Fig. (2]-(b), we show S* obtained by the two methods mentioned at the beginning of 
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FIG. 2: (a): Temperature dependence of So- (b): Temperature dependence of S* obtained from 
real frequencies (filled circles) and from Matsubara frequencies (open circles). Square lattice. 



Sec. IV The solid circles represents S* by fitting ImL n (uj) and ImL 12 (w) in real frequency 



at tu — > oo limit. The open circles represent S* computed using Eq. (|27[) and Eq. (28). The 



values of S* at open and closed circles are very close, indicating the consistency between the 
real and Matsubara frequency approach to calculate S*. 

The dependence on electron density of So an d S* is more non-trivial, which is difficult to 
tell from analytical formulas. Fig. [3] shows Sq and S* at various densities for U = 1.75D. Sq 
changes sign from positive at half filling to negative as electron density decreases, while S* 
remains negative. The behavior of So here is also due to the breakdown of coherence as the 
evolution of spectral weight. In a doped Mott insulator, the quasiparticle peak gradually 
diminishes as the system is doped away from half-filling [27]. Thus near half-filling, the 
transport is dominated by the coherent excitations near Fermi surface. But when the doping 
is heavy enough to kill quasiparticles, transport is carried by incoherent excitations in the 
Hubbard bands. Therefore So turns to a same sign with S*, since S* is dominated by the 
Hubbard bands(see Fig. [l}(c) and discussion there). In Fig. [3] we also put the results of S* 
by real and Matsubara frequency approach. 



B. Triangular lattice 

Recent interest on thermoelectric performance of correlated systems was attributed to 
the discovery of TEP enhancement in highly electron doped cobaltates|28j. The Co atoms 
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U=1.75D, T=0.125D 
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1- 




0.75 0.8 0.85 0.9 
electron density 



0.95 



FIG. 3: Doping dependence of So and S* for U = 1.75D. S* was obtained from real frequencies 
(filled circles) and Matsubara frequencies (open circles). The temperature here is T = 0.125D. 
Square lattice. 

in the CoO-i layers form a triangular lattice. The physics behind the large TEP in Na x Co02 
is highly non-trivial. For example, the Na potential is crucial to induce the correlation in 
Nao.7CoC>2$9\, and the spin and orbital degrees of freedom are argued to be a key factor 
for the enhancement [51 [30]. These complexities are beyond a single band Hubbard on 
a triangular lattice. Here we only focus on some qualitative features of Sq and S* in a 
electron-doped single band Hubbard model on triangular lattice. 

In triangular lattice, U = 12\t\ and we use a positive t. S* in this section is solely 



computed by using Eq. (27) and Eq. (27). 



Fig.|3j-(a) and -(b) shows the density dependence of So and S* for two different interaction 
strength. Here we present the full range for electron doping. Here S* is from the summation 
over Matsubara frequency. For U = 1.25D(Fig. [41(a)), So is negative near half- filling and 
changes to positive after a small amount of doping. As the density approaches to band 
insulator(ra = 2), the merging of Sq and S* is very evident. For smaller interaction strength, 
i.e., U = 0.5-D, 5*0 and S* also display similar trend through the range of electron density. 
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FIG. 4: (a) and (b): Density dependence of S and S* for U = 1.25D and J7 = 0.5D. S* was 
obtained from the Matsubara frequencies. Triangular lattice. 

This behavior is similar to the case on square lattice, Fig. [3} The discrepancy between So 
and S* is most evident for U = 1.25D and near half-filling (n = 1.0), since around this 
regime the coherent quasiparticles take a significant role in transport. For electron density 
larger than 1.5, which is the range of interest for cobaltate, the trend of S* shows that it is 
a reasonable approximation to Sq. 



V. SUMMARY 

Using the formulae derived in Sec. [TTJ we investigate to what extent the AC limit of 
thermoelectric power, S*, can be a reasonable approximation to the DC limit, Sq. Analytical 
and numerical results on a single-band Hubbard model show that below and around coherent 
temperature, i.e., when the spectral weight around quasiparticle peak dominates in the 
thermoelectric transport, the behaviors of So and S* are significantly different. Specifically, 
Sq displays multiple sign changes around the coherent temperature, but S* does not. But 
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when the temperature is well beyond the coherent regime, thus the transport properties are 
dominated by the incoherent excitations, S* shows same sign and similar magnitude to So 
and can give reasonable prediction on the behavior of S*. 



Our work suggest that a realistic implementation of Eq. (27) and Eq. (28 ) in LDA+DMFT 
codes can serve as a useful guide for the search of high performance thermoelectric materials 
among the strongly correlated electron systems, which have a very broad temperature regime 
characterized by incoherent transport. 

At the time of writing, we are aware of a recent work by M. Uchida et al.[3I], in which the 
incoherent thermoelectric transport over a wide temperature range is studied in a typical 
density-driven Mott transition system La\- x Sr x VO^ and the validity of Mott-Heikes formula 
for real strongly correlated materials is verified. 
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Appendix A: Sum rules for ReLi2(w) and ReLn(w) in DMFT 

In this appendix, we compute L\ x and L* l2 m the framework of dynamical mean field 



theory and show they also obey the general formulae, Eq. (31 and Eq. (32) 
In terms of retarded current-current correlations, 



ReL7?(w) = --Im 



(Al) 



which can be computed in Matsubara frequencies by standard diagrammatic techniques [ 
In the infinite dimension limit, a significant simplification is achieved because all nonlo- 
cal irreducible vertex collapse and only the first bubble diagram survives p21 [21]. This 
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simplification leads to 



2 



f(co')-f(co' + co) 



A k (u')A k (u' + to). 



(A2) 



Now we calculate L*-. Using Eq. (15) 



Changing variables by 

LO\ = OJ + J 

leads to 

2 



x ( fW-fi» + ^ '\ MunMu/ + u) . (A3) 



to 



fe,<7 ^ Xj J 

+ 2 EffrO [ duxdub-^fMAuiuJAkfa). (A4) 



The sum rule J dcoiA k (ui) = 1 simplifies the first term to 



In the second term, Kramer-Kronig relation can be used to eliminate the integral over u\, 
i.e., 



Then we use the fact that 



and 



/ du 1 ^^- = -ReG k (u 2 ). 
J 



2ReG k (u)lmG k (uj) = ImG?(w) 
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to simplify the second term on the right hand side of Eq. ( A4 ) to 



Applying integration by part over k, it turns out to be 



J du 2 u 2 f{u 2 )A k {u 2 ). 



de k \ 2 , / d 2 e k 



Combined with the first term, we have 

The calculation for L\ 2 is similar and straightforward, which results in 



f(u)A k (u). (A6) 
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